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ABSTRACT 

The first detection of tfie silicate absorption feature in AGNs was made at 9.7 
for the prototypical Seyfert 2 galaxy NGC 1068 over 30 years ago, indicating the 
presence of a large column of silicate dust in the line-of-sight to the nucleus. It is now 
well recognized that type 2 AGNs exhibit prominent silicate absorption bands, while 
the silicate bands of type 1 AGNs appear in emission. More recently, using the Mid- 
Infrared Interferometric Instrument on the Very Large Telescope Interferometer, Jaffe 
et al. (2004) by the first time spatially resolved the parsec-sized dust torus around 
NGC 1068 and found that the 10 /im silicate absorption feature of the innermost hot 
component exhibits an anomalous profile differing from that of the interstellar medium 
and that of common olivine- type silicate dust. While they ascribed the anomalous 
absorption profile to gehlenite (Ca2Al2Si07, a calcium aluminum silicate species), we 
propose a physical dust model and argue that, although the presence of gehlenite is 
not ruled out, the anomalous absorption feature mainly arises from silicon carbide. 
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1 INTRODUCTION 

Amorphous silicate is a major dust component of the Galac- 
tic interstellar medium (ISM), as revealed directly by the 
smooth, featureless 9.7 /xm Si-O stretching and 18 /im O- 
Si-0 bending vibrational absorption bands, and indirectly 
by the depletion of the silicate-forming elements Si, Mg, Fe 
and O in the gas phase (see Draine 2003). The first detec- 
tion of the silicate absorption feature in active galactic nuclei 
(AGNs) was made at 9.7 ^m for the prototypical Seyfert 2 
galaxy NGC 1068 35 years ago by Rieke & Low (1975) and 
Kleinmann et al. (1976), indicating the presence of a large 
column of silicate dust in the line-of-sight to the nucleus. It 
is known now that most of the type 2 AGNs (which only 
exhibit narrow emission lines in their optical spectra) dis- 
play silicate absorption bands, while the silicate features of 
type 1 AGNs (which are characterized by the presence of 
both broad and narrow emission lines in their optical spec- 
tra) appear in emission. This is consistent with the unified 
theory of AGNs which suggests that essentially all AGNs 
were "born equal": as illustrated in Figure la, all types of 
AGNs are surrounded by an optically thick dust torus and 
are basically the same object but viewed from different lines 
of sight (see e.g. Antonucci 1993; Urry & Padovani 1995). 

Within the framework of the unified theory of AGNs, 
for type 1 AGNs which are viewed face-on, one would ex- 
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pect to see the silicate features in emission since the silicate 
dust in the surface of the inner torus wall will be heated to 
temperatures of several hundred kelvin to ~ 1000 K by the 
radiation from the central engine, allowing for a direct detec- 
tion of the 9.7 ^m and 18 /im silicate bands emitted from the 
dust (see Figure la). Their detection has been reported with 
Spitzer (Hao et al. 2005; Siebenmorgen et al. 2005; Sturm et 
al. 2005; Weedman et al. 2005; Shi et al. 2006; Schweitzer et 
al. 2008). For type 2 AGNs which are viewed edge-on, the 
silicate features are mostly seen in absorption (see Figure 
la) because of the obscuration of the optically thick dust 
torus (e.g. see Roche et al. 1991, 2007; Siebenmorgen et al. 
2004; Hao et al. 2007; Spoon et al. 2007). 

We note that the silicate emission features have also 
been detected in several type 2 QSOs (Sturm et al. 2006, 
Teplitz et al. 2006). Efstathiou (2006), Marshall et al. (2007) 
and Schweitzer et al. (2008) suggested that the silicate emis- 
sion may actually arise from the narrow-line region (NLR) 
dust. However, Nikutta et al. (2009) argued that the clumpy 
dust torus model could explain the detection of the sili- 
cate feature in emission in type 2 sources. More recently, 
by combining Gemmi and Spitzer mid-IR imaging and spec- 
troscopy of NGC 2110 (the closest known Seyfert 2 galaxy 
with silicate emission features). Mason et al. (2009) were 
able to constrain the location of the silicate-emitting region 
to within 32 pc of the nucleus. Their results were consistent 
with both a NLR origin and an edge-on torus origin of the 
silicate emission. 
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Figure 1. (a): A schematic illustration of the dust torus around AGNs. The unified model of AGNs predicts the silicate feature to 
appear in emission in type 1 AGNs (which are viewed face-on), and in absorption in type 2 AGNs (which are viewed edge-on), (b): A 
schematic illustration of the innermost ~ 2pc dust structure of the near nuclear region of the prototypical Seyfert galaxy NGC 1068. The 
interferometric observations of Jaffe et al. (2004) resolved this dust structure into two components: a well-resolved 2.1pcx3.4pc warm 
component and an inner hot component of ~0.7pc along the jet and < 1 pc across, (c): A comparison of the Galactic center silicate 
absorption optical depth profile (Kemper et al. 2004) with the interferometric mid-IR spectrum of the innermost ~ 2 pc near nuclear 
dust region of NGC 1068 (Jaffe et al. 2004). 



The silicate emission profiles of AGNs are rather di- 
verse, suggesting that the AGN silicate grains were proba- 
bly not "born equal". Some AGNs display an "anomalous" 
Si-0 stretching band which is considerably shifted to longer 
wavelengths (from the canonical 9.7 /im to ~ 11.5 /xm; here- 
after we call it "redshift" ) and broadened in comparison with 
the Galactic silicate profile (Hao et al. 2005, Siebenmorgen 
et al. 2005, Sturm et al. 2005, Smith et al. 2010), with the 
degree of redshifting and broadening varying from one AGN 
to another. 

Some AGNs, e.g. 3C273 (Hao et al. 2005) and 
NGC 7213 (Wu et al. 2009), exhibit an unusually strong red 
tail of the 18 /^m 0-Si-O bending band, incomparable with 
the Galactic silicate profiles. The silicate absorption pro- 
files of some AGNs also show appreciable deviations from 
that of the Galactic ISM. More specifically, using the Mid- 
Infrared Interferometric Instrument (MIDI) on the ESO's 
Very Large Telescope Interferometer (VLTI), Jaffe et al. 
(2004) by the first time spatially resolved the parsec-sized 
dust torus around NGC 1068 (see Figure lb) and found that 
the 10 /^m silicate absorption feature of the innermost hot 
component exhibits an anomalous profile differing from that 
of the ISM (see Figure Ic): the 9.7 /xm silicate absorption 
spectrum of NGC 1068 shows a relatively fiat profile from 
8 to 9 /im and then a sharp drop between 9 and 10 /im, 
while the Galactic silicate absorption profiles begin to drop 
already at about 8 /im. Jaffe et al. (2004) found that the 
profile of the silicate absorption towards the hot compo- 
nent does not fit well to the profiles of common olivine-type 
silicate dust (also see Raban et al. 2009). They obtained a 
much better fit using the profile of gehlenite (Ca2 Al2Si07), a 
calcium aluminum silicate species (a high-temperature dust 
species found in some supergiant stars, e.g. see Speck et al. 
2000). 

As a part of a systematic study of the nature of the 
dust in AGNs and particularly of the AGN dust mineralogy, 
we quantitatively model the anomalous silicate absorption 
feature associated with the innermost hot component of the 
dust torus of NGC 1068. 



2 NGC 1068 

At a distance of only 14.4 Mpc, NGC 1068 is one of the clos- 
est and probably the most intensely studied Seyfert 2 galaxy. 
It has played a key role in the establishment of the unified 
model of AGNs: classified as a Seyfert 2 based on the pres- 
ence of narrow emission lines and absence of broad emission 
lines in its optical spectrum, NGC 1068 is known to also har- 
bor an obscured Seyfert 1 nucleus as revealed by the detec- 
tion of broad emission lines in polarized light (Antonucci & 
Miller 1985). As expected from an edge-on geometry of the 
proposed dust torus for type 2 AGNs, NGC 1068 has long 
been known to display a strong absorption band at 9.7 /im, 
attributed to amorphous silicate dust (Rieke & Low 1975, 
Kleinmann et al. 1976, Roche et al. 1991, Lutz et al. 2000, 
Le Floc'h et al. 2001, Siebenmorgen et al. 2004). 

Jaffe et al. (2004) obtained the 8-13.5 /im A*'-band spec- 
tra of the nuclear region of NGC 1068, using the MIDI in- 
strument coupled to the VLTI [with a spatial resolution 
of ~ 10 milliarcsec (mas) at A=10/im, corresponding to 
~0.7pc]. Their interferometric mid-IR observations by the 
first time spatially resolved the central parsec-sized circum- 
nuclear dust structure of NGC 1068 (see Figure lb). 

Rhee & Larkin (2006) obtained spatially-resolved mid- 
IR spectra of the nucleus of NGC 1068, using the Long 
Wavelength Spectrograph (7.71-12.48 /im) of the Keck I 10 m 
telescope (with a spatial resolution of 0.25", correspond- 
ing to ~18pc in physical scale). They found that the sil- 
icate absorption feature varies over the nuclear region: the 
depth of the feature reaches its maximum in the nucleus and 
decreases with distance from the central engine. Spatially- 
resolved mid-IR spectra of the nucleus of NGC 1068 were 
also obtained by Mason et al. (2006), using the Michelle 
A''-band (7-14 /im) spectrometer of the Gemini North 8.1m 
telescope (with a spatial resolution of 0.4", corresponding 
to ~ 30 pc). They also found that the silicate feature profile 
and depth exhibit striking spatial variations. 

The interferometric observations of Jaffe et al. (2004) 
resolved the near-nuclear mid-IR emission into a warm com- 
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Figure 2. Comparison of the interferometric spectrum of the 
innermost, parsec-sized dusty nuclear region of NGC 1068 (Jaffe 
et al. 2004) with the model spectra calculated from various dust 
species: (a) SiC (red line) or amorphous olivine silicate dust (green 
line) for both the inner hot component and the surrounding warm 
component; (b) gehlenite in the hot component and SiC in the 
warm component (red line), or SiC in the hot component and 
gehlenite in the warm component (green line); (c) amorphous 
olivine in the hot component and SiC in the warm component 
(red line), or SiC in the hot component and amorphous olivine 
in the warm component (green line); and (d) gehlenite in the hot 
component and amorphous olivine in the warm component (red 
line), or gehlenite in both the hot and warm components (green 
line). 



ponent (~ 320 K) in a (well-resolved) structure ~ 2.1 pc thick 
and ~ 3.4 pc in diameter, surrounding a smaller (marginally- 
resolved) hot structure (> 800 K) of ~0.7pc in size along 
the jet and < Ipc across (also see Raban et al. 2009). They 
found that the silicate absorption feature could not be fit 
with the Galactic silicate profile or the profiles of common 
olivine- type silicate dust (see Figure Ic), as they begin to 
drop already at ~ 8 fim, while the interferometric spectrum 
of NGC 1068 (with the highest spatial resolution) shows a 
relatively fiat profile from 8 to 9 fim (and then a sharp drop 
between 9 and 10 /im). Therefore, they concluded that their 
observations would require special dust properties in the in- 
nermost 2pc around the nucleus of NGC 1068. They ob- 
tained a quite satisfactory fit with the profile of gehlenite 
Ca2Al2Si07, a high-temperature dust species. However, the 
approach taken by Jaffe et al. (2004) seems to us too simpli- 
fied (see §5). This stimulates us to present a more physical 
model (see §3) to interpret the unusual silicate absorption 
spectrum of NGC 1068 obtained by JaflFe et al. (2004). 

3 THE MODEL 

Following Jaffe et al. (2004), we assume that the nuclear 
dust region of NGC 1068 consists of two components: a hot 
component and a warm component. Let and be 
the hydrogen column densities of the warm and hot com- 
ponents, respectively; and Th be the (mean) dust tem- 
peratures of the warm and hot components, respectively; 
C'abs('^i'^) and C^ba(fli-^) be the absorption cross sections 
of the dust species of size a at wavelength A of the warm 
and hot components, respectively. For simplicity, we will as- 
sume a single grains size of a = 0.1 /im for both the warm 
and hot components. Let N^^^^ and A'^^^t be the column 



densities of the dust (of size a) of the warm and hot com- 
ponents, respectively. Let o"^s(A) and o-^bs(-^) be the total 
absorption cross sections per H nucleon of the dust at wave- 
length A of the warm and hot components, respectively. Ap- 
parently, we have a^y,,{X) = C^b,ia,X) x (N^^JN^) and 
'^abs(A) = Ci,(a,A) X (K^^JN^). Let and fi^ be the 
solid angles extended by the warm and hot components, 
respectively. The flux emitted by the warm and hot compo- 
nents received at the Earth is given bjQ 

^"-B,{T^) {1 - exp [-N^a:Um 



Fx = 



47J- • ■ - - - - 

+ ^BA(r„) {l - exp [-Ka^U^)] } 
xexp[-7VHO-™bs(^)] • 



(1) 



The exp [— TVg (j!|"bg(A)] term in eq.2 accounts for the fact 
that the emission from the inner hot component is further 
absorbed by the dust in the warm component (like a screen). 
The solid angles of the warm and hot components are taken 
to be Qw = 30 mas x 49 mas ~ 3.45 x 10"^'' sr, and Qh ~ 
10 mas X 12 mas x 2.82 x 10"^^ sr (Jaffe et aL 2004). 

We take the column density of the dust of a given 
species to be limited by the cosmic abundance of the rarest 
element contained in a given species (e.g. Ca in gehlen- 
ite Ca2Al2Si07). Therefore, the total number of gehlenite 
grains (per H) can not exceed 



gch 



[Ca/H] 



Ah 



(2) 



2 (47r/3) a^pgeh/ (^gchmn) 

where [Ca/H] is the Ca abundance (relative to H), pgah = 
2.91 gcm~^ is the mass density of gehlenite (Mutschke et 
al. 1998), /igch = 274 is the molecular weight of gehlenite 
(Ca2 AbSiO?) , and mn is the atomic mass of H. Similarly, for 
amorphous olivine-type (MgFeSi04) silicate dust, the total 
numbers of silicate and SiC grains (per H) are limited to 

A,ii _ [Si/H] A,i, _ [Si/H] 

Ah (47r/3) a^p^jj/ (^^jj^h) ' Ah (47r/3) a^paic/ (/iaicmn) 

where [Si/H] is the Si abundance (relative to H), p^w = 
3.5gcm~^ is the mass density of MgFeSi04, psii ~ 172 is 
the molecular weight of MgFeSi04, psic ~ 3.2 gcm"'^ is the 
mass density of silicon carbide dust, and /Xsic = 40 is the 
molecular weight of SiC. 

In modeling the interferometric mid-IR absorption spec- 
trum, we first need to specify the dust compositions for the 
warm and hot components. Once the dust compositions are 
specified, we obtain their absorption cross sections from Mie 
theory (Bohren & Huffman 1983). We are then left with only 
four parameters: T™, Th, A^h , and A^. 

4 RESULTS 

We have considered a wide variety of dust materials, includ- 
ing amorphous olivine (Dorschner et al. 1995), amorphous 



^ Consider a dusty system with a hydrogen column density of 
A'^H- Let o"^i3s(A) be the total dust absorption cross section per H 
nucleon. Let px be the power radiated per solid angle in [A, A + dA] 
per H. With the dust self-absorption taken into account, the 
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Figure 3. Comparison of the interferometric spectrum of tlie in- 
nermost, parsec-sized dusty nuclear region of NGC 1068 with the 
model spectra calculated from (a) a mixture of 80% amorphous 
olivine silicate dust and 20% amorphous carbon; (b) a mixture of 
80% SiC and 20% amorphous carbon; (c) a mixture of 50% SiC 
and 50% amorphous olivine; and (d) a mixture of 90% SiC and 
10% gehlenite. We assume the same dust composition for both 
the hot and warm regions. 



pyroxene (Jager et al. 1994), amorphous carbon (Rouleau & 
Martin 1991), a-SiC (Laor & Draine 1993), glassy gehlen- 
ite (Mutschke et al. 1998), and spinel (Fabian et al. 2001). 
We have also considered mixtures of these materials, with 
their absorption properties calculated from Mie theory in 
combination with the Bruggeman effective medium theory 
(Bohren & Huffman 1983). For the abundance constraints 
(see eqs.3,4), we adopt the solar abundance of [Ca/H] ~ 
2.2-10-'^ and [Si/H] ^ 31.6-10"^ (Asplundet al. 2009). With 
a = 0.1/xm, iVgci/A'^H ~ 4.1-10~", iYsii/Afn ~ 6.2-10"", and 
iVsic/A''H ~ 1.6 ■ 10-" (see eqs.2,3). 

As shown in Figure 2, the best fits to the interfero- 
metric mid-IR spectrum of Jaffe et al. (2004) are provided 
by models which all require SiC (see Table 1 for the model 
parameters): either with pure SiC in both the hot region 
and the surrounding warm region (see Figure 2a), or with 
gehlenite in the hot region and SiC in the warm region (see 
Figure 2b) , or with amorphous olivine in the hot region and 
SiC in the warm region (see Figure 2c). The composition of 
the inner hot component is less well constrained as the large 
column density required for the inner component causes the 
dust there to emit effectively like a blackbody. The absorp- 
tion profile is predominantly provided by the SiC dust in 
the outer warm region. 

We have also tried pure amorphous olivine dust for both 
regions. As shown in Figure 2a, even the best fit fails to fit 
the observed spectrum, particularly in the blue wing. With 
the warm region filled with amorphous olivine and with the 
hot region filled with either SiC (see Figure 2c) or gehlenite 
(see Figure 2d), we obtain model spectra similar to that 
of the pure amorphous olivine model, confirming that it is 
the dust in the outer warm region dominates the absorption 
profile. 

We have also considered gehlenite. Assuming either 
gehlenite for both regions (see Figure 2d) or gehlenite only 
for the outer warm region and other dust species (e.g. SiC) 
for the inner hot region (see Figure 2b), we see that good 
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Table 1. Model Parameters. The dust column densities A^dust can 
be obtained from A^h using eqs.2,3 (e.g. N^^i = A''h x [Afgd/Afn] )■ 
The dust masses m^ust can be derived from A^dust a-nd the solid 
angle Q, (e.g. rrigei =0.5 [Ca/H] A'gei^H.^^H x f^i^^)- 



fits to the blue wing of the 9.7 ^m feature are achieved, but 
both models fail at A > 11 ^m. 

Finally, we model the observed mid-IR spectrum with 
a mixture of two dust species and assume the same mix- 
ture for both the hot and warm regions (see Figure 3). For 
the amorphous olivine-amorphous carbon mixture, even the 
closest fit (given by a mixture of silicate of a volume frac- 
tion of 80% and amorphous carbon of a volume fraction of 
20%) fails to reproduce the blue wing of the 9.7 /im feature 
(see Figure 3a). Better fits could be achieved if the mixture 
contains SiC (see Figure 3b,c,d). 

5 DISCUSSION 

In §4 we have shown that, for a physically plausible model, 
to achieve any reasonably close fits to the interferometric 
spectrum of NGC 1068, SiC is required. SiC is a highly re- 
fractory dust species. It has been identified in the outfiows of 
carbon stars. Interstellar SiC grains have been found intact 
in "primitive " meteorites (e.g. Bernatowicz et al. 1987), al- 
though the fraction of silicon in the local diffuse ISM in SiC 
dust is small: < 5% (Whittet et al. 1990) or ~ 3% (Min et 
al. 2007). SiC was one of the three major dust species (i.e. 
amorphous silicate, graphite, and SiC) considered by Laor 

6 Draine (1993) in their extensive AGN dust modeling. Al- 
though we are not able to rule out gehlenite as the major 
dust species in the inner hot component, SiC should be the 
dominant dust species in the outer warm component and 
dominates the mid-IR absorption. 

For all models, the determined temperatures of the hot 
and warm components are > 800 K and < 300 K, respec- 
tively (see Table 1). For the SiC-containing model, the tem- 
peratures of the dust in the hot and warm components are 
Th ~ HOOK and ~ 250 K, respectively. The derived H 
column density of the hot component A^h is in the order of 
~ 10^^ cm"'^, consistent with the measurements of the inner 
wall by Matt et al. (1997, 2004). The warm component re- 
quires to be in the order of ~ 10^^ cm"^, consistent with 
the measurements of Kinkhabwala et al. (2002). 

Finally, we note that although the calculations were car- 
ried out for a grain radius of 0.1 /im, increasing the grain 
radius or assuming a porous structure causes the deviation 
between the observed and model spectra to increase. 

Jaffe et al. (2004) argued that a quite satisfactory fit 
to the interferometric mid-IR spectrum of the near-nuclear 
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dust region of NGC 1068 could be obtained with the absorp- 
tion profile of gehlenite. Gehlenite is a silicate of the melilite 
group with a high condensation temperatures (~ 1735 K for 
melilite, see Hutchison 2004). It has been found in objects in 
the solar system: calcium- aluminum rich inclusions in car- 
bonaceous chondrites (see e.g. Morlok et al. 2008), anhy- 
drous interplanetary dust particles (Greshake et al. 1996), 
and cometary dust collected by the Stardust spacecraft from 
comet Wild 2 (Zolensky et al. 2006, Schmitz et al. 2008). It 
is assumed that gehlenite was produced close to the young 
Sun at high temperatures and was distributed in the entire 
system. 

In their model, Jaffe et al. (2004) assumed that (1) both 
the warm and hot components emit like a black-body, and 
(2) each component is subject to an attenuation of exp(— ta), 
with the optical depth profile t\ varying to match that of 
known silicate species. Their modeling theme can be de- 
scribed as the following 

Fx = c„Ba(T») exp(-rr)+ChBA(rh) exp(-r^) exp(-rr) , (4) 
where Cm and Ch are the constants (related to the total dust 
masses) for the warm and hot components, respectively^ 
Tu, and Th are the dust temperatures of the warm and hot 
components, respectively; Bx{T) is the Planck function of 
temperature T at wavelength A; rjj" and are the optical 
depths in front of the warm and hot components, respec- 
tively. Since the hot component is embedded in the warm 
component, the hot component is subject to an attenuation 

While it was rather successful in fitting the MIDI spec- 
tra, the assumptions made by Jaffe et al. (2004) were not 
justified: (1) their modeling theme implied that as if there 
were an absorbing screen in front of the hot component and 
an additional absorbing screen in front of the warm compo- 
nent; and (2) if the grains in the hot and warm components 
indeed emit like a blackbody as they assumed, the grains 
must be much larger than several micrometers and reach 
the geometrical optics regime0 These grains would not pro- 
duce any noticeable spectral features around A = 10 fim, 
and therefore the silicate absorption must arise from the 
absorbing screens. 

To summarize, we agree with Jaffe et al. (2004) that the 
interferometric spectrum of NGC 1068 cannot be explained 
with amorphous silicate with an olivine normative composi- 
tion. However, unlike Jaffe et al. (2004) who attributed the 
unusual 9.7 ^m absorption feature of NGC 1068 to gehlen- 
ite, we argue that SiC is mostly responsible for this feature 
although the presence of gehlenite in the inner hot region is 
not ruled out. 
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